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ABSTRACT 

It is well established that supermassive black holes in nearby elliptical galaxies correlate tightly 
with the kinematic property (Mbh— cr correlation) and stellar mass (Mbh — Mhost correlation) of their 
host spheroids. However, it is not clear what the relations would be at the low-mass end, and how 
they evolve. Here, we investigate these relations in low-mass systems (Mbh ~ 10 6 — 10 8 M Q ) using 
the Aquila Simulation, a high-resolution cosmological hydrodynamic simulation which follows the 
formation and evolution of stars and black holes in a Milky Way-size galaxy and its substructures. 
We find a number of interesting results on the origin and evolution of the scaling relations in these 
systems: (1) there is a strong redshift evolution in the Mbh - f relation, but a much weaker one in 
the Mbh~ Mhost relation; (2) there is a close link between the Mbh— o relation and the dynamical 
state of the system - the galaxies that fall on the observed correlation appear to have reached virial 
equilibrium. (3) the star formation and black hole growth are self-regulated in galaxies - the ratio 
between black hole accretion rate and star formation rate remains nearly constant in a wide redshift 
span z — — 6. These findings suggest that the observed correlations have different origins: the 
Mbh— o relation may be the result of virial equilibrium, while the Mbh— Mhost relation may the result 
of self-regulated star formation and black hole growth in galaxies. 

Subject headings: scaling relations- low mass black holes - host galaxies - fundamental plane - 
coevolution methods: numerical - hydrodynamical 



1. INTRODUCTION 

A major development in observational astrophysics in 
recent years is the discovery that most, if not all, nearby 
elliptical galaxi es host a supermassive black ho le (SMBH) 
at their center ( Kormendy fc Richstone| [ l99 5;), and that 
the masses of the SMBHs correlate tightly with the global 
properties of the spheroid components of their hosts, 
such as the s tellar velocity dispersion (the M B h~ o cor- 
relation, e.g., |Ferrarese fcHVIcrritt 2000; Gcbhardtetal 



20001 |Tremaine et al.||2002| |Gultekin et al.||2009b| [ Gra- 
ham et al.| 2011 ) , and the stellar masses (the Mbh M SP h 



corre l ation, e.g., [Magorrian et al.|1998||Marconi fc Hunt| 



2003 Haring & Rix 2004[ ). These correlations suggest 
that the formation and evolution of SMBHs and their 
host galaxie s are closely linked (e.g., Haehnelt & Kauff- 
mann|2000|). 

'ML J 



The origin of these scaling relations, however, remains 
a hot debate. A number of different models have been 
proposed to explain the correlations, such as feedback 

Silk fc Rees| 1998l 



from active galactic nucl ei (AGN, e.g 



Ciotti & Ostriker 2007), gas co mpetition between star 
formation and BH accretion (e.g., Li et al.|2007al|Jahnke| 
fc Macci6|2011 ), galaxy mergers (e.g., Feng|2UU7 [^and a 
combinat ion of merger and AGN feedback ( e.g., Li 



2007b[ [Hopkins et a .J[2006a| [20071 |2009a|b[ ). 

In order to unravel the origin of the correlati 



et al. 



tions, recent 

observational efforts have focused on their evolution in 
hope to identify the crucial physical processes i n galaxy 
formation which give rise to the relations (e.g., [Shields 
20031 |Treu et al.||2004| |Walter et al.||2004| [S hields 
'20061 IMcLure et al.||2006[ |Peng et al.||2006[ [ Woo 
alviander et al. 120071 ITreu et al. 120071 fWoo 



et al 
et al. 



et al. 



2006 



qxzl25@psu.edu ~ 



et1QTl[2008l |Jahnke et aLl[20 09; Bcnn ert et al.||2010| [De^] 
carli et al.|2010||Merloni et al.|2010||Cisternas et al.|2011 [ 



Bennert et al. 2011). These surveys suggest that there 
is either weak or no evolution in the Mbh— & relation up 
to z ~ 2, but there is a strong redshift evolution in the 
Mbh— M sp h relation. However, it was also suggested that 
the relations between Mbh and total host-galaxy lumi- 
nosity or stellar mass, Mbh— Mhost, may not be evolving, 
or at least not as rapidly as the Mbh— M B ph relation (e.g., 
Jahnke et aTl|2009l JBennert et all |2lH0} [Merloni et al. 
2010| |Cisternas et al.|[2011[ |Bennert et al.||2011[ ). 

These results suggest that the Mbh— M sp h and Mbh— o 
relations have different origins. More importantly, the 
difference between Mbh— M sp h and Mbh— Mhost provides 
a crucial clue on the origin and evolution of the mass 
scaling relation, because M B h— M sp h is closely connected 
to the formation of b ulges, which is believed to result 



from major mergers ([Barnes fc Hernquist 



quist 



1992 



1992 
and 



thus 



Hcrn- 

He"- 



j ,19931 |Hopkins et al.| 2006a| ), 
pends str ongry~on redshift and t rie merging history of 
a galaxy (Hopkins et al. 2010a|b[ ). On the other hand, 
M B h— Mhost may simply reflect the growth of stars and 
BHs in a galaxy regardless its type. Therefore, the 
Mbh - Mhost may represent a more fundamental prop- 
erty of galaxies. In this work, we adopt the latter as it 
is also easy to measure from the simulations. 

Theoretically, the evolution of the BH - host rela- 
tions has been studied using numerical simul ations (e.g., 
Robertson et al . 2006 ; Di Mattco ct al. 2008 ; |Johansson| 
et al.||2009t I Booth fc Sc haye 2011 D and semi-analytical 




However, these models differ m their predictions lor the 



evolution of the scaling relations. In particular, |Robert-| 
son et al. ( 2006 ) employed a set of idealized merger simu- 
lations at various redshifts and found that the slope of the 
Mbh - f relation remains roughly constant at redshifts 
z = — 6, but the normalization shows a weak redshift 



dependence. |Di Matteo et al. ( 2008 1 performed the first 
direct cosmological hydrodynamic simulation with BHs 
down to redshift z — 1, and found a weak redshift evolu- 
tion in the normalizatio n and slope of both correlations. 
Booth & Schaye ( 2011[ ) also used a cosmological simula- 
tion but a different BH model and found that, while they 
reproduced the local correlations and the observed evolu- 
tion of Mbh— Mhost for massive galaxies, they predicted 
an evolution in the Mbh - a which contradicts with cur- 
rent observations. 

These studies have painted a confusing picture of the 
BH - host correlations. The confusion may come from 
the different samples used, and the different evolution- 
ary stages of these galaxies. The galaxy samples used 
to derive the original correlations are dominated by local 
massive ellipticals and inactive galaxies. Recently, it has 
been suggested that the correlations va ry with galaxy 
type and mass. Gultekin et al. (2009b) showed differ- 
ent scaling relations tor ellipticals and for spirals, and 
the latter have a larger int rinsic scatter in the Mb h~Q" 



relation than the ellipticals. Kormendy et al. (2011 ) sug- 
gested that the BH masses do not correlate with galaxy 
pseudobuldges or disks. Moreover, |Lauer et al.| p007a[ ) 
found that, at the very high mass end, the BH masses 
derived from Mbh— g differ those from Mbh— M sp h by 
nearly one order of magnitude. 

At the low mass end (Mbh < 1O 8 M ), most of the ob- 
served galaxies are active galactic nuclei (AGN). It was 
reported that, compared to inactive galaxies, AGNs show 
a larger dispersion in the Mbh— c relation, and that the 



lipticals (e.g., Greene et al.||2008 


2010 


Kuo et al. 112011 


Bennert et al. 201 1| Xiao et al. 


201 1|. More interest- 



a large deviation from the "classical" Mbh— o relation, 
they seem to f ollow the Mbh— M sp h relation. Recently, 



Decarli et al. ( |2012[ ) compiled a sample of 26 quasars 
in the low-mass end (Mbh ~ 10 7 — 1O 9 M ), and found 
that the Mbh— M sp h relation is consistent with that of 
the quiescent galaxies. 

In order to understand the origin of the scaling rela- 
tions, it is important to follow the formation and evo- 
lution of BHs and their host galaxies, in particular at 
the low-mass end, because these are the building blocks 
of the massive ones at the present day. Based on our 
current understanding of structure formation in the cold 
dark matter (CDM) cosmology, the central BHs may fol- 
low the same hierarchical assembly of the galaxies as well. 
Black holes in the mass range of 10 5-8 M Q thus serve as 
a bridge between the SMBHs and their seeds, thus would 
provide crucial information on the growth of BHs. 

In this work, we study the BH - host relations and their 
evolution in low mass systems using a high-resolution 
cosmological hydrodynamic simulation which includes 
important physics of dark matter, gas dynamics, star 
formation, black hole growth, and feedback pro cesses. 
The simulation used the Aquila initial condition ( |Scan-| 
napieco et al. 20121, which was constructed to produce 



a Milky Way-size halo and its sub-structures (referred to 
as "Aquila Simulation" hereafter). We also explore the 
origin of the relations by comparing model predictions 
with observations. The paper is organized as follows: in 
§2, we describe our models and computational methods. 
We present the evolution of the relations in §3, models 
for their origins in §4, and summarize in §5. 

2. THE AQUILA SIMULATION 

In order to achieve high resolutions in a cosmologi- 
cal simulation to study the BH - galaxy correlations at 
different cosmic time, we carried out the Aquila Simu- 
lation, which follows the formation and evolution of a 



Milky Way-size g a laxy and its substructures (Wadepuhl 
& Springel 2011 Scannapieco et al. 2012 1. The ini- 



tial condition is orig inally from the Aquarius Project 



( Springel et al.| 2008 ) , which produced the largest ever 
particle simulation of a Milky Way-sized dark matter 
halo. The hydrodynamical initial condition is recon- 
structed from the original collisionless one by splitting 
each original particle into a dark matter and gas particle 
pair fWadepuhl fc Springel||2011[ ). 

The Aquila Simulation includes dark matter, gas dy- 
namics, star formation, black hole growth, and feedback 
processes. Star formation is modeled in a multi-phase 
ISM, with a rate that follows the Schmidt-Kennicutt Law 
( Schmidt||1959 Kennicutt 1998). Feedback from super- 
novae includes both thermal and kinetic forms. Thermal 
feedback is captured through a multi-phase model of the 
ISM by an effective equation of state for star-forming 
gas, and the k inetic feedback is modeled a s a galactic 
wind based on Springel & Hernquist (2003). We adopt 



a constant wind velocity of t> w ind = 484 km s , a mass- 
loss rate that is twice of the star formation rate, and 
an energy efficiency of unity such that the wind carries 
100% of the supernova energy. The wind direction is 
anisotropical, preferentially perpendicular to the galactic 
disk. This wind model causes an outflow of gas, trans- 
porting energy, matter and metals out of the galactic disk 
in proportion to the star formation rate. 
The model of black h o le grow th an d feedback follows 



that o f Springel et al. (2005a) and Di Matteo et al 



(20051, where black holes are represented by collision 



less sink" particles that interact gravitationally with 
other components and accrete gas from their surround- 
ings. The accretion rate is estimated from the local gas 
density and sound speed using a spherical Bondi-Hoyle 
( Bondi|[l952) iBondi'fc Hoyle||1944j |Hoyle fc Lyttleton 
1941[ ) model that is limited by the Fddington rate. Feed- 
back from black hole accretion is modeled as thermal 
energy, ~ 5% of the radiation, injected into surrounding 
gas isotropically. This feedback scheme self-regulates the 
growth of the black hole and has been demonstrated to 
successfully reproduce many observed properties of local 
elliptical galaxies (e.g,, [S pringel et al. 2005a|b Di Matteo 
et al. 2005 Hopkins et al. 2006a]) and the most distant 
quasars at z ~ 6 ( |Li et al.||2007b[ ) . Mergers of black 
holes can happen once two black holes are close enough 
and their relative speed is less than the local sound speed. 



that of previous work (Di Matteo et al.|2008 


Sijacki et al. 


2009 


Di Matteo et al. 


2012): a seed black 


hole of mass 



iBH = iu- it "1V10 

tial minimum of each new halo identified by the friends- 
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Figure 1. Evolution of the BHs and galaxies from the Aquila Simulation. The images show the projected density of both gas (color-coded 
by temperature: blue indicates cold gas, brown indicates hot, tenuous gas) and stars (represented by the bright yellow color). The black 
holes are represented by the red dots, the size of which is proportional to the BH mass. Feedback from both supernovae and accreting BHs 
creates hot bubbles around the galaxy centers. The box size is 10 h~ 1 Wlpc in comoving coordinates. 



of-friends (FOF) group rinding algorithm with a total 
mass greater than 10 10 /i -1 M©. 

The computation was performed using the parallel, N- 
body/SPH code GADG ET-3, which is an imp r oved ver 



sion o f that described in Springel et al. (2001); Springel 
(2005). For the computation of gravi tational fo rces 



thecode uses the "TreePM" method (|Xu||1995 ) that 
combines a "tree" algorithm (Barnes & Hut 1986 1 for 
short-ran ge forces and a Fourier t ransform particle-mesh 



method ( jHockney &: Eastwood 
GADGET i mplements 



forces. 



1981 1 for long-range 



he entropy-cons erving 



formulation of SPH ( |Springel fc Hernquist] |2002[) with 
adapti ve particle smoothing, as in Hernquist &; Katz 



(1989). Radiative cooling and heating processes ar e cal- 



culated assuming collisional i onization equilibrium ( Katz 
et al.||1996[ |Dav<§ et al.||1999D, and the UV background 



model of Haardt & Madau (1996) is used, which describes 
a spatial uniform UV background leading to reionization 
roughly at z ss 6 in the simulation. 

The whole simulation falls in a periodic box of 
100 /i~ 1 Mpc on each side with a zoom-in region of a size 
5x5x5 h~ 3 Mpc 3 . The spatial resolution is ~ 250 hr 1 
pc in the zoom-in region. The mass resolution of this 
zoom-in region is 1.97 x 10 5 /i _1 M for dark matter par- 
ticles, ~ 1.875 x 10 4 fe _1 MQ for gas and star particles. 
The cosmological parameters used in the simulation are 
Q m = 0.25, Qa = 0.75, cr 8 = 0.9 and h = 0.73, consistent 
with the five-year results of the WMAP (Komatsu et al. 



2009). The simulation evolves from z = 127 to z = 0. 

In the simulation, each snapshot is processed by an on- 
flying FOF algorithm which tests the dark matter linking 
length as if they are less than 20% of their mean spac- 
ing. Gas and star particles are then linked to the nearest 
dark matter particle. A substructur e detection algorit hm 
SUBFIND (an extended version of IDolag et al.||2009| for 
gas and star particle) is then applied to each group by cal- 
culating the local density and searching for locally over- 
dense region. Such an overdense region will be marked 
as a substructure of the parent group. Throughout this 
work, a galaxy is defined as the group returned by the 
SUBFIND, which includes dark matter halo, gas, stars, 
and black holes. We only select galaxies from the high- 
resolution, zoom-in region. 

Figure [l] shows the evolution of the BHs and galaxies 
from z ~ 6.2 to z = from the Aquila Simulation. The 
gas follows the distribution of dark matter in filamentary 
structures, and stars form in high density regions along 
the filaments. The most massive galaxy in the zoom- 
in region resides in the intersection of the filaments, the 
highest density peak in the simulated volume where gas 
concentrates in the deep potential well. The BHs form in 
these massive halos, and they grow through gas accretion 
and mergers following the hierarchical buildup of their 
host galaxies. The most massive BH resides in the main 
halo at the present day, and there are about a dozen of 
BHs in the mass range of 10 5 — 10 8 M Q . 
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Figure 2. The Mbh~ ct relation at different rcdshifts from the 
Aquila Simulation, in comparison with the local observation. The 
red dots represent the simulation data, the red solid curve is the 
fitting of our sample using a power-law formula, while the black 
solid c urve is the best fit of the nearby galaxies by Trcmainc ct al. 
(2002), and the dashed lines indicate the range of the fitting. 



3. THE BH - HOST RELATIONS AND THEIR EVOLUTIONS 

3.1. The Mbh~ct Correlation 

From the Aquila Simulation, we select a sample of 
galaxies with a stellar mass of > 10 8 M Q , which con- 
tains over 5000 star particles. This criterion enables ro- 
bust measurement of the galaxy properties. Most of the 
selected galaxies are dwarfs, and have actively accreting 
BHs. For each galaxy, the BH mass is computed di- 
rectly in the simulation, and we calculate the projected 
half-mass effective radius R e , and the line-of-sight mass- 
wieghted stellar velocity dispersion a within R e . This 
projection procedure is performed for over 100 random 
line-of-sights. 

The resulting Mbh - o~ relation is shown in Figure [2] 
at six representative redshifts, in comparison with the 
observed M ^—a rel ation of the nearby galaxies by 



Tremaine et al. (2002 1. We fit our data at each redshift 



with t he sa me power-law formula as that in Tremaine] 



et al. (20021: 



log 



M 

~Mr. 



BH 



alog( 



200 kms- 



where a and j3 represent the slope and the normalization 
of the relation, respectively. The best fit of the local 
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Figure 3. The evolution of the parameters of the Mbh - g rela- 
tion, slope a (in black) and normalization (in red), with redshift 
in the simulation. The solid lines are least-squares fittings of the 
data. 



Mbh— c relation from |Tremaine et al. ( 2002 1 gives a slope 
of a = 4.02±0.32, and a normalizationoFF= 8.13±0.06. 

From Figure [2] there is a clear trend of evolution of 
the Mbh - cr relation in both normalization and slope. 
The modeled galaxies systematically lie below the local 
Mbh~ f relation at high redshift. A correlation emerges 
quickly after the BHs are formed, as high as z — 4, which 
is also the fast growth phase of BHs. As the galaxies 
evolve to low redshift, they converge to the observed lo- 
cal relation quickly. By 2 = 0, some dwarfs eventually 
reach the line defined by local massive ellipticals, but 
the overall slope is noticeable steeper. The steeper slope 
at the low-mass end was also seen in simulations of IDil 
Matteo et al. ( 2008 ) , and in low-mass AGN observations 



oT Greene et al. ( |2010 ) from precise BH mass measure 



2011 



ments using megamaser (note, however, Xiao et al. 
reported a shallower slope with a large sample of aGm s. ) . 
The change in the Mbh~ ct relation at low mass proba- 
bly indicates a different mode of BH growth in these ob- 
jects c ompared with more m assive galaxies with classical 



bulges Greene et al. (2008). 



The evolution of the Mbh— a relation, in form of slope 
a and normalization j3, is shown in Figure [3j The scat- 
tering in both slope and normalization becomes larger 
at higher redshift, due to smaller galaxy number and in- 
creasing difficulty in the measurement of a. Nevertheless, 
both slope and normalization show strong evolution with 
redshift out to z — 6, as illustrated by the fitting lines: 



a = 4.38 - 0.212, = 8.07 - 0.20z 



(2) 



The evolution of the normalization of the Mbh~ g re- 
lation from our simulation is consistent with the result 



of Di Matteo et ah] ( 2008 ) from cosmological s imulations, 
which also s howed cx — 0.20z, and that of |Robertson| 
et al. (2006) from merger simulations, oc — 0.186z by 



(1) fixing the slope t o 4.0, b ut shallower than th at given by 



Booth fc Sc haye 
et al.| ( [2006) aM 



42011b, 
L)i Mat 



cx — 0.32z. B oth Robertson 



Matteo et ah (2008) suggested that 



the expected velocity dispersion for a given stellar mass 



4 



is larg er at higher redshift. A s for the evolution of the 



slope, Robertson et al 
trend possibly due to t 



(20061 did not see a monotonic 
le limited number of redshifts in 
that study. Instead they put an upper limit o n the abso- 
lute value o f the evolution of slope, |</>| < 0.3. Di Matteo 
et al. ( |2008| ) showed that the slope becomes steeper with 
decreasing redshift, similar to the trend we find, except 
that they have a much larger slope at z ~ 3 — 4 since the 
massive BHs (> 10 8 M ) fall above the local correlation. 
The steep "tilt" at z ~ 3 — 4 in their work may be caused 
by the rapid growth of these massive BHs which gained 
their masses rapidly through mild super-Eddington ac- 
cretion rates during this period. In fact, a similar trend 
is also present in Figure [3] at z ~ 3 — 1 when the most 
massive black holes in our simulation assembled most of 
their masses. 

In observation, the Mbh-ct relation has been studied 
by using the width of some emission lines suc h as QUI 
or H a as a proxy for stellar v elocity dispersion ( Nelson & 



Whittle 
this sea 



1996||Gu et al.|2009 l. These studies sugg est that 
mg relation either does not evolve (e.g., [Shields 



et al. 2003), or does so weakly, with BHs more massive 
than interred fro m the local relation by a factor of a few 



at z 



1 (e.g., Woo et al. 2008 Gu et al. 20091. In 



Figure [3j the evolution ot the Mbh <t relation at z < 1 
is much weaker than that indicated by the above fittings. 
At z — 0, our resul t is within the observed range of low - 
mass BH systems ( Greene et al.|[2010 Xiao et aT~||2011 ). 



3.2. The M B H-M host Correlation 

For the same galaxy sample as in the previous section, 
we measure the masses of the hosts. Currently, we do 
not have an effective procedure to decompose the galaxy 
from a cosmological simulation into bulge and disk com- 
ponents, thus we are using half of the galaxy stellar mass 
as a proxy of bulge mass. The resulting Mbh— Mh os t re- 
lation is shown in Figure |4j in comparison with obser- 
vations. We fit our data at each redshift with the same 



power-law formula as in Haring &: Rix ( 2004 ) 



, M BH , , 
l0g M^ al0g( " 



Af, 



10 n M 







-) + b 



(3) 



where a and b represent the slope and the normalization 
of the relation, respective ly. The best f i t of t he local 
MBH-Mhost relation from Haring & Rix (2004) gives a 



slope of a — 1.12 ± 0.06 and a normalization of b 
8.20 ±0.10. 

Compared to Mbh— f relation in Figure [2| the 
Mbh— Mhost relation shows a weaker evolution with red- 
shift. Most of the BHs in our simulation have masses 
Mbh < 10 7 M Q , and they lie below the local correlation, 
while the more massive ones Mbh > 1O 7 M , lies above 
the line, which in broa d agreement with simulations of 



Di Matteo et al. (20081, and observations which have a 



above and below the local relation I 


e -fr? 


Greene et al. 


2010 


Bennert et al.||2011 


Xiao et al. 


2011 





The evolution of the IV 



iBH - Mhost relation, m rorm ot 
slope a and normalization 6, is shown in Figure [5j The 
scattering in both parameters becomes larger at nigher 
redshift, due to smaller galaxy count, but on average, 
both slope and normalization do not show strong evo- 
lution with redshift out to z — 6, as illustrated by the 
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Figure 4. The Mbh - Mhost relation at different redshifts from 
the Aquila Simulation, in comparison with the local observations. 
The red dots represent the simulation data, the red solid curve is 
the fitting of our sample using a power-law formula, while the black 
solid c urve is the best fit of the nearby galaxies by |Haring &: Rix| 
| |2004| l, and the dashed lines indicate the range of their httmg. 
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Figure 5. The evolution of the parameters of the Mbh— Mhost 
relation, slope a (in black) and normalization b (in red), with red- 
shift in the simulation. The solid lines are least-squares fittings of 
the data. 



5 



fitting lines. If we just focus on z < 2, there is no notice- 
able evolution in either the slope or the normalization. 



1.38- 0.04z,6 = 7.64 + 0.02z 



(4) 



The modest evolution in the Mbh— Mh os t relat ion from 
our simulation a grees well with the result from |Di Mat-' 



teo et al. (|2008j) in the same mass range. It appears 
that the BHs in massive galaxies (M * > 10 10 Mp) grow 



more quickly than the stell ar mass (Di Matteo et al. 
2i)(),S| IBooth fc Schaye||201 1j^ This may - be due to the 



super-Fddington accretion ( Di Matteo et al. 1 20081) and 
the increased effici ency early on (Booth & Schayc; {2011| ) 



in these models. Johansson et aT~(2009p argued that 



SMBHs are less likely to develop before their parent 
bulges as they found that over massive BHs were not 
evolving to the local correlations f rom t heir merger sim- 
ulations. However, Naab et al. (2009) suggested that 
minor mergers could be a solution to the increase of stel- 
lar masses. In our simulation, most of the BHs are low 
mass, and they accrete at sub-Eddington rates most of 
the time, which differs significantly from the massive ones 
in the previous studies. 



tion has been studied by a number of groups (e.g., [Shields 


et al. 2003; Treu et al.|2004 


Peng et al. 12006; Woo et al. 


2006 


Salviander et al.|2007 


Treu et al.||2007[|Woo et al. 


2008 


IJahnke et al. 20091 i 


iennert et al. 20101 Decarli 


et al.|2010; Merloni et al.|201 


0||Cisternas et ai.|2011||Ben- 


nert et al.|2011). In particu 


ar, it was recently suggested 



masses of the host spheroids M sp h evolves stro ngly with 
M BH /M sph oc (1 + z) 196 (|Bennert 



time out to z 



et al. 2011). However, most of the studies derived the 
bulge stellar mass from the luminosity of the spheroid. 
While this can be done robustly for local ellipticals, the 
uncertainty becomes larger at higher redshift, because 
the stellar populations in bulges are not well understood 
at z > 1, it is more difficult to distinguish between the 
host and the bulge component, and bulge may not have 
formed in some high redshift galaxies or AGNs. In fact, 
when considering the total luminosity or the stellar mass 
of the host galaxy, it was found that the Mbh— Mh os t has 
a mu ch weaker evolution , if any, than the Mbh~ M sr ,h 
(e.g., IJahnke et al.|[2 009l |Bennert et al.||2010[ |Merloni 
et al.||2010[ |Cistern as 
consistent with our 
Studies of quasar 



asHjTalTpTllj |Bennert et al.||201l| ) 
finding in Figure [5] 
galaxies 



that 



host galaxies indicate 
Mbh— Mhost rela tion may be present as early as z w 4 
( Peng et al.|2006 l, and that evidence for evolution of this 
relation is not conclusive because of the following three 
reasons: (1) BH mass calibration has a factor of 2 system- 
atic uncertainty; (2) bul ge mass is not directl y measured; 
and (3) selection bias. 



(2010) studied 96 



Decarli et al. 

quasars and their host galaxies and found an evolution 
in Mbh— Mg a i from z = 3 to the present day. However, 
the possibility that these host galaxies may be biased to- 
wards lower m ass galaxies at high r edshift is present, as 



pointed out by( Lauer et al 
lead to a more massive blac. 



2007b[ ). Such a bias would 
c hole for a give luminosity 
by nearly one order of magnitude. On the other hand, 
a limi ted sample of obscured AGNs from |Sarria et al.| 
(2010) suggests that Mbh— Mh os t relation still holds at 
high redshift. 
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Figure 6. Ratio of BH accretion rate to star formation rate of 
our galaxy sample at different redshift from the Aquila Simulation. 
The black error bar indicates the range of the ratio at a given time, 
while the black dot represents the mean v alue. The red and blue 
dots represent observational data based on iGultekin et all l|2009a|l 
and |Hainline et al,| l |2012| |, respectively. The red, dasned line is the 
least-squares htting to the simulation data. 

The difference between M sp h and Mhost provides an 
important clue on the origin and evolution of the scal- 
ing relations, because the local Mbh— o and Mbh— M sp h 
are determined by the formation of bulges, which are be 



lieved from ma 


or mergers (e.g., 


Hopkins et al. 


2006a). Overall, 


Mbh - f reflects the change of 



TJTe 



[1992 

evolution 



1993 
of the 



while M B h— Mhost simply reflects the growth history of 
a galaxy and its central BH, regardless its type and or 
activity. 

4. ORIGIN OF THE SCALING RELATIONS 

The origin of the BH scaling relations has been an open 
question. Our simulation shows that there is a remark- 
able evolution with redshift in the Mbh— c relation, but 
little or no evolution in the M B h— Mhost relation. The 
different evolutions of the scaling relations suggest that 
they may have different origins. 

4.1. Mbh - Mhost Relation as a Result of Self-regulated 
Growth in Galaxies 

In order to understand the mass correlation, we turn 
to the growth history of the galaxies in our sample. Our 
simulation includes a number of physical processes of 
baryonic matter, including gasdynamics, cooling, star 
formation, BH accretion, and feedback from both SNe 
and BHs. Dark matter halo quickly forms and signifi- 
cant amount of mass is gathered through merger, lead- 
ing to the reconfiguration of the gravitational potential 
and relaxation of DM particles. Gas falls into forming 
gravitational potential well and is shocked to high tem- 
perature. A considerable fraction of gas cools, loses its 
angular momentum, pile up in the center of the halo and 
stellar component builds up a short period of time, as 
well as the central existing BH. Feedback from SNe and 
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AGNs heats up gas to hot phase, compensate the cold gas 
inflow and regulate both star formation and BH growth. 
The BHs grow until self-regulated, when its feedback en- 
ergy can sufficiently unbind in-falling gas and halt accre- 
tion as well as star formation. The difference between an 
isolated merger and a cosmological simulation is that a 
constant gas supply from IGM is present in a cosmolog- 
ical simulation, which will exhibit as a cooling inflow to 
the galaxy in the early stage. 

As a result of the interplay among these physical pro- 
cesses, the star formation and BH growth in galaxies ap- 
pear to be self-regulated and balanced. Figure [6] shows 
the ratio of BH accretion rate to star formation rate of 
each galaxy at different redshift in the simulation, in 
comparison with two available observational data points 
at z = and z ~ 2, respectively. The lo cal galaxy sam- 
ple is b ased on X-ray lumino sity from jGultekin et al.| 



(2009a), H a luminosity from |Ho et al.| (|2003p and in- 



frared luminosity from Melendez et al.| (|2u(J8|) . The SFR 



is calculated from H a luminosity or fa r infrared luminos 

(fl998l). 
ith 



ity using the conversion relations by Kcnnicutt 



The accretion rate is estimated with the same method as 

We convert X-ray luminosity 



Mullaney et al. (2012) 



to bolometric luminosity with a constant factor of 22.4, 
and assume a radiation efficiency of 0.1 to calculate the 
accretion rate. The total number of the galaxies i n this 



sample is 12. Th e high redshift AGN sample from |Hain 
line et al. (2012) includes 11 galaxies with both SFR and 



AGJN bolometric luminosity. We calculate the accretion 
rate using the same method above. 

As shown in Figure [6j although there is a large dis- 
persion in the data due to variation in the BH accre- 
tion rate and star formation rate, the least-squares fit- 
ting to the simulation data (the red, dashed line) shows 
log(M BH /SFR) = -2.754 - 0.02z. The fitting indicates 
nearly negligible evolution in the ratio over redshifts 
z = — 6, and it lies close to the available observational 
data. More interestingly, this growth ratio, Mbh/SFR. = 
—2.754, is very close to the mass ratio of the observed 
M RH -M hnHt relatio n, M BH /Mhost ~ (1 - 2) x lO" 3 (e.g., 
Hating fc Rix1[2004| . 

This result suggests that BHs have grown in-step with 
their host galaxies since at least z ~ 6, that the BH 
growth and star formation in galaxies are self-regulated, 
probably owing to feedback and availability of gas sup- 
ply, and that the ratio of the two remains constant re- 
gardless the triggering mechanism. The fact that the 
ratio maintains constant at Mbh/SFR = —2.754 natu- 
rally gives rise to the observed Mbh - Mh os t relation in 
the local universe. 



Recently, Mullaney et al. (2012) studied the growth of 
SMBH and stars in a large sample of galaxies in the red- 
shift range 0.5 < z < 2.5. They found that the average 
SMBH accretion rate follows remarkably similar trend 
with stellar mass and redshift as the average SFR of 
their host galaxies, and that the ratio Mbh/SFR ~ 1 0~ 3 
similar to o ur finding in this work. Furthermore, |Shi| 
et al. ( 2009| ) studied infrared properties of 57 SDSS type 
-1 quasars at z~l and found that the average ratio of 
Mbh/SFR in quasar hosts shows little evolution with 
redshift. These observations give strong support to our 
hypothesis that the observed M B h — M host relation is a 
result of self-regulated and balanced growth in galaxies. 
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Figure 7. The ratio of kinetic energy to the potential energy of 
stellar components within half-mass radius of the galaxy at differ- 
ent redshift from the Aquila Simulation. This ratio serves as an 
indicator of the dynamical state of a galaxy, for example, if the 
ratio is 1, it means the system is in virial equilibrium. The black 
error bar indicates the ratio range at each redshift, the black dot 
represents the mean value, the red dot is the value derived from 
observations of Marconi & Hunt (2003|, the yellow dot in dicates 
data f rom Sani et al. ( 2U11 ), and blue d ot indicates d ata from|Jiang| 
|et al.|pQll[ |. Tne samples of |Marconi fe Hunt| | |2003^ and |Sam et al.| 
| |2U11^ are dominated by massive elliptical galaxies which sa tisfy 
the local Mbh-" correlation, while that of Jiang et al. (20fl) are 
low-mass systems with BHs of 10 5 — 10 6 Mq , which show significant 
deviation from the local Mbh - f correlation. 

Our model predicts that the Mbh - Mh os t relation is a 
global, fundamental property of galaxies, it does not de- 
pend on galaxy type or mass,, and it does not evolve with 
redshift. However, it does not rule out strong evolution 
in the relation with spheroid, because the Mbh - M sp h de- 
pends on the formation of bulges, which is believed to be 
triggered by major mergers of galaxies, therefore, proper- 
ties of spheroids depend sensitively on environment and 
redshift which determine the merging history of a galaxy. 

4.2. Mbh - c Relation as a Result of Virial Equilibrium 

We have shown in the previous section that the galaxies 
from our simulation fall on the the observed local Mbh - c 
relation at z = 0, but there is a remarkable evolution in 
the relation with redshift, while observations of massive 
galaxies with spheroids show little evolution. In order to 
understand the origin of this relation and the cause of 
the evolution and discrepancy, we turn to the dynamical 
state of the galaxies as it determines the stellar velocity 
dispersion. 

A number of model s have been propos ed to explain 



the Mbh - o relation. Silk & Reesl (fl998| assumed the 



formation of SMBH proceeds galaxies and the energetic 
wind from quasars expels the gas, resulting in a regula- 
tion of the growth of th e spheroidal stellar component. 
Ciotti & Ostriker (2007) argued that radiative feedback 
from BHs in elliptic als give rise to this scaling relation. 
Hopkins et al. ( 2007 ) suggested that it takes a dynamical 
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time for the feedback energy from the BHs to unbind the 
surrounding gas and reached a fundamental-plane, relat- 
ing the BH mass with a combined parameter with stellar 
mass and velocity dispersion. 

The finding that the Mbh— 17 and the fundamental- 
plane relation are closely related suggest that galaxies 
satisfying these relation s may have reached relaxation, or 



gested by Binney & Tremaine (2008 1 : 



virial equilib rium (e.g., Robertson et al. 2006 Hopkins 
et al.|2009a l, which satisfies the virial theorem: 



2K + W = 



(5) 



where K = 3/2M st cr 2 is the kinetic energy of the stars, 
and W = — 3/5GMM*/Re is gravitational potential of 
the galaxy with mass M within the effective radius R e . 

To examine the dynamical state of our galaxies, we de- 
fine a parameter virial ratio, A = 2K/W. If the system 
is in virial equilibrium, then A = 1. Figure [7] shows the 
resulting virial ratio A of the galaxy sample at different 
redshift from our simulation, in comparison with avail - 
able observational data from Marco ni fe Hunt (2003), 



Sani et al. 



_ J2°11 Land | Jiang et ai | ( |2011 - 
Sani et al.| ( |201l| ) studied the black hole-host galaxy 
scaling relations m Spitzer/IRAC band. The velocity 
dispersion, K-band luminosity are taken from their Ta- 
ble 2 and the effective radius R e from their Table 3. The 
galaxy mass is then derived from the K-band luminos- 
ity with a mass-to-light ratio M/L = 1. The number 
of total gala xies from this sample is 50. In |Marconi~fe" 
Hunt ( |2003[ ), since the effective radius R e was measured 
in J band, we use the J band luminosity to compute 
stellar mass also assuming M/L = 1. This gave us to- 
tal 37 gala xies. Note the sa mples of |Marconi fc Hunt] 
(2003) and Sani et al. (2011) are dominated by massive 
elliptical galaxies whi ch satisfy th e local Mbh 



(2011) 



a cor- 
are low-mass 



relation, while that of Jiang et al 
systems with BHs of 10° — 10° M Q , which fall below the 
local Mbh~ o correlation. A sample of galaxies with low 
mass BHs is also included as we cross check the ones with 
both photometric measurements from Jiang et al. fl2011|) 
and v elocity dispersion measurements from |.X.iao et al.| 
(2011). Luminosity of the bulge is estimated from the 
surface brightness in I band a nd the luminosity d istance 

The 
This 



(2011) 



inferred from redshift z from Jiang et al 
ratio M/L in this band is also assumed to be 
gives us a sample of 30 galaxies. 

From Figure [7J a clear evolution of the virial ratio A is 
seen: the mean value flattens around ~ 1.1 — 1.2 up to 
z ~ 2, then increases gradually with redshift and reaches 
~ 2.2 at z ~ 6. At z = 0, our galaxies agree with the ob- 
served Mbh— o relation, and our mean value of the ratio 



coni & Hunt 


(2003 


) and 


Sani et al. 


(2011 



This agree- 



-a correlation are indeed in virial 



mcnt is exciting, beef 
that fall on the Mbh~ 
equilibrium . On t he other hand, most of the galaxies of 
Jiang et al. (2011 1 fall below the classical Mbh— a corre- 
lation, like ours at z > 3, therefore, it is not surprising 
that its 2K/W ratio is larger than 1, similar to our result 
at z ~ 6. This means that the galaxies at high redshift 
in our simulation have not yet reached virial equilibrium. 

This result suggests that the Mbh— o correlation may 
be a result of virial equilibrium of the galaxies. Quali- 
tatively, from the criterion for virial equilibrium as sug- 



< o A >~ 0.45- 



GM 



R r 



(6) 



which gives a relation between the galaxy mass and the 
stellar velocity dispersion and effective radius: 



M oc a z R P 



(7) 



Since the luminosity of a galaxy, L (oc M*), is propor- 
tional to the stellar surface brightness and the projected 
surface area: 

L cx R 2 J e (8) 

Assuming a constant mass-to-light ratio M/L and a 
constant surface brightness / e , the above two relations 
lead to 

L(cx M.) cx a 4 (9) 

From our results in the previous section, the self- 
regulated growth of stars and BHs in galaxies gives 
Mbh oc M*, therefore, 



M BH oc M* 



(10) 



This virial equilibrim-origin model can explain the dis- 
crepancy between observations and our simulation. In 
the observations, most of the galaxies, in particular mas- 
sive ones and those with bulges, may have reached relax- 
ation and virial equilibrium quickly, because the more 
mas sive a galaxy is, the shor ter relaxation timescale it 
has ( Binney & Tremaine 2008 ) . It is believed that merger 
can trigger a rapid growth of B Hs and the build up o f 
ellipticals and bulge component (Springel et al. 2005a), 
when large amount of shocked gas by gravitational in- 
teraction fuels BHs accretion and intense star formation 



(Springel et al. 2005a ). This scenario is als o supported 
by ob s ervati ons (e.g., Bennert et al. ( 2008 ) ; Alexander 
et al. (2005)). The spheroids formed through violent 
merger quickly relax and reach virial equilibrium. As 
a result, they satisfy the Mbh— o relation and show little 
or no evolution. On the other hand, low-mass systems 
take a longer time to reach virial equilibrium, thus they 
may show some evolution in the Mbh— c- 

5. DISCUSSIONS 

In order to achieve desired resolutions to resolve low- 
mass BH systems in cosmological simulations, we focus 
on a Milky Way-size halo in this paper. However, this 
has led to a serious limitation: the galaxy sample is very 
small. The large scattering due to the small number 
statistics in our results, in particular at high redshift, 
may have affected the actual evolution of the scaling re- 
lations we studied. 

In such a zoom-in simulation, there might be con- 
tamination by heavier low-resolutio n particles outside o f 
the high-resolution zoom-in region ( Springel et al.|2008[ ). 
This could be more severe as we approach to the bound- 
ary of the zoom-in region. However, most of our galax- 
ies are located well within the boundary. We checked 
the presence of low-resolution particles in all the subha- 
los which host BHs, and found that the mass fraction is 
very small, about 10 -3 — 10~ 4 . We conclude that the 
contamination is not significant in our galaxy sample. 
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The seeds of galactic BHs is an outstanding problem. 
It is unknown where the seed of the Milky Way came 
from. In our s imulation, we follow the previous studies Li 



et al.| (2007b); Pi Matteo et al.| ( |2008| and create seeds of 
10° M Q in halos with mass above 10 M Q . Although it is 
not our focus to reproduce exactly the same BH as in the 
Milky Way, our seeding scheme produced a ~ 10 s M BH 
in the main halo, much more massive than the current 
estimate of 10 6 M Q in the Milky Way. 

A wide range of s eed masses, from 10 M to 1O 6 M , 
has been proposed (Volonteri 2010). For examples, a 
10 2 Mf?> by the collapsed remnants from 



small seed of 



the first stars (Abel et al. 2002; Bromm & Larson 2004 


Gao et al. 2007 


Yoshida et al.|2008|), or a massive one o 



~ 10" M by the cat astrophic collapse of a supermassive 
star, or gas clump (ICarr et aTT] | 1984| i Bromm k, Loeb 



2003} IBegelman et al. 2006; Tanaka & Haiman 20IM7" 

Recently, | Volonteri fc INatarajan| ( |2009[ ) studied the ef- 
fects of BH seeds on the BH scaling relations using semi- 
analytical merger trees. They investigated two seeding 
models, light seeds of 10 2 M© from PopIII stars, and 
heavy seeds from direct gas collapse (which has a distri- 
bution function with a peak of 1O 5 M ), and suggested 
that the Mbh— c relation is a result of seeding mecha- 
nism and growth prescription, and that massive seeds 
produce better agreement with observations at z = 0. 
They also predicted a large population of low-mass BHs 
at high redshifts, which has not been detected. 

We plan to improve the current study in future work 
by increasing the statistics. We will perform uniform 
cosmological hydrodynamic simulations in larger boxes 
in order to obtain a large number of galaxies at different 
redshift. We also plan to investigate different BH seeding 
schemes and the resulting BH - host galaxy relations. 

6. SUMMARY 

We have investigated the Mbh— cr and relations 
Mbh~ Mhost relations in low-mass systems (Mbh ~ 
10 6 — 1O 8 M ) using the Aquila Simulation, a high- 
resolution cosmological hydrodynamic simulation fo- 
cused on a Milky Way-size galaxy. The simulation in- 
cluded dark matter, gas dynamics, star formation, black 
hole growth, and feedback from both SNe and accreting 
BHs. Here is a list of our findings: 

• The Mbh - cr and Mbh - Mhost relations evolve dif- 
ferently with redshift: the former shows a remark- 
able evolution in both slope and normalization, 
while the latter shows little evolution in the red- 
shift range z = — 6. 

• There is a close link between the Mbh— cr relation 
and the dynamical state of the system - the galax- 
ies that fall on the observed correlation have a virial 
ratio A = 2K/W close to 1, indicating that they are 
in virial equilibrium, while those that off the rela- 
tion instead have a ratio larger than 1. 

• The star formation and black hole growth in galax- 
ies are self-regulated - the ratio between black 
hole accretion rate and star formation rate remains 
nearly constant over z = — 6, probably owing to 
feedback and gas availability in the galaxies. 



These results suggest that observed scaling correlations 
have different origins: the Mbh— cr relation may be the 
result of virial equilibrium, while the Mbh— Mhost rela- 
tion may the result of self-regulated star formation and 
black hole growth in galaxies. 
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